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1 Introduction 

Let VL be an open bounded set in R^, with a Lipschitz boundary. We consider the 
following turbulent circulation model : 

-div(z/(A;)VM) = / in Vt 
-(l\Y{a{k)Vk) = v{k)\Vu\^ m Vt 
M = on dVt 
k = on d^l 



(P) 



Here /, a and u are given, and the functions u,k : Q ^ are the unknowns. 
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We study Problem (P) under the following main assumption: 



/ G with r > I 



< a, z/ : M"*^ — s> M"*" are continuous 
3 5>0: a{s),i^{s) >5 Vs G M 



+ 



Problem (P) is a simplified scalar version of the RANS model arising in oceanography 
(see [10,11,1]): the function u is an idealisation of the mean velocity of the fiuid and 
k is the turbulent kinetic energy. The mathematical analysis of (P) is a step towards 
better understanding the RANS model. Various studies were made in this direction. 
Some existence results were established in [10,7]. 

In this paper we focus on the case where the viscosity functions a and u are not a 
priori bounded. In fact (see [11,7]), in the relevant physical situation, we have 



We will establish an existence result for a weak solution for (P) under less restrictive 
assumptions than in [7]. An important feature is that our assumptions are satisfied 
under (Hp), contrarily to the assumptions made in [7]. 

Moreover we give additional regularity results for the weak solution we obtain. In 
particular, under (Hq) and the following additional assumption: a is proportional to 
u, dQ is of class C^'", / G C^'"{Q) and u G C^'"(]R~^), we prove the existence of a 
classical solution for (P). 

We also compare our results with the results presented in [10]. 

Another feature of our work is to considere various notions of weak solution for 
Problem (P): H^-solution, if-solution, distributional solution, renormalized solution, 
'energy solution', classical solution. We give some relations between these notions. 

1.1 Notions of weak solution for (P) 

We can reformulate equation (P).2 by using the Kirchoff transform. Let 





Instead of (P).2, we can consider 



(P).2' 



= voA-^{K)\Vu\'^ onfi 



where K = A{k). 

In fact, from every distributional solution K G of (P).2' we obtain a distri- 

butional solution k of (P).2 by setting k = A~^{K). This property is related to the 
facts that A is invertible, A~^{0) = and |y4~^(s)| < Cs (this can be seen by using 
the assumptions made on u in (Hq)). 

The situation is more complicated for equation (P).l, where the a priori unbounded 
coefficient z/(/c) appears in the principal part of the operator and cannot be removed. 
Hence we have to restrict m to satisfy the energy condition 

/ z/(A;)|Vn|2 < oo. (1) 

Nevertheless we will see later on that various non equivalent notions of weak solution 
can be considered for (P).l. 

We will introduce the notions of W-solution and H-solution. It is also possible to 
consider the notion of renormalized solution (see [10] chap. 5). In [7] the authors 
defined another notion that they call energy solution. 
We will give some relations between these notions in the Appendix I. 

Remark now that under the restriction (1), the right hand side in (P).2 (or in 
(P).2') is only a priori in L^{Q). Hence (see [2]) it is natural to seek k in the space 

np<3/2W^0''(^^)- 

We want to find a function u vanishing on dQ that satisfies the energy condition 
(1). This leads to considering the following spaces: 

Wk = !^veH'oi^)-h']k<^] 

Hk = closure of C^{^) with respect to [.]k 
where we used the notation 

For any measurable function k, the map defines a norm on Wk- In the general 
situation and Wk are not equal. Moreover Wk is not necessarily complete and 
a function in does not always have a uniquely defined gradient (see [14]). If we 
assume that z/(A;) G L^{Q) then Wk is complete and in fact C Wk are Hilbert 
spaces (see [6,14,13]) when they are equipped with the scalar product 

{v,w) = / i'{k)Vv\/w. 
Jn 



Consequently, we will consider the following two distinct notions of solution for 



(P).l: 



u is called a iffe-solution of (P).l if u G Hk and 
^ u{k)VuVv = fv Wv e Hk 



u is called a Wfc-solution of (P).l if m G Wk and 



Finally, we define the following notions of weak solution for (P): 
(u, k) is called a H-solution of (P) if 

k e np<3/2Vro'^(l^), u e Hk, 

k is a distributional solution of (P).2 and n is a i/fc-solution of (P).l 
{u, k) is called a W-soIution of (P) if 

k e np<3/2<''(^), u^Wk 

is a distributional solution of (P).2 and m is a VF^-solution of (P).l 
1.2 Main results 

Let {Hi) and {H2) denote the following conditions: 

{Hi) 37 > : a{s) > -fu{s) Vs G M+ 
{H2) 37 > : a{s) = 7i/(s) Vs G M+ 

We will establish: 

Theorem 1 Assume that (Hq) and (Hi) hold. Then there exists at least one W- 
solution {u,k) for (P) such that 



Then the W-solution {u,k) given in Theorem 1 is a distributional solution of (P). 

Theorem 2 Assume that (Hq) and (H2) hold. Then the W-solution {u,k) given in 
Theorem 1 satisfies 




(2) 



Corollary 1 Assume that in addition to (Hq) and (Hi) we have 



3 1^0 > : z/(s) < z/o(l + s^), Vs G M+. 



(3) 



G C°'°(l]), /or some a G (0,1). 



(4) 



Moreover [u, k) is also a H-solution of (P) (and in fact a classical weak solution). 
If in addition to (Hq) and (H2) we assume that dVL is of class C^'", / G C°'"(f2) and 
V e C^'"(M+) then 

u,keC^'^{U), for some p e (0,1), 
and {u, k) is a classical solution of (P). 



1.3 Discussion of the results 



In Theorem 1 we give an existence result of a W-solution. We next give some regu- 
larity results: firstly the property (2) and secondly (in Theorem 2) the property (4). 
Finally, in Theorem 2 we give an existence result for a classical solution for (P). 



The main previous studies of Problem (P) are presented in [10] chap. 5 and in [7]. 



In [10] chap. 5, the authors prove the existence of a renormalized solution for (P) un- 
der the assumptions (i^o) and (-^2)- It seems that their proof also works under {Hq) 
and {Hi). Nevertheless the notion of renormalized solution is very weak. A renor- 
malized solution (u, k) for (P) is a distributional solution if z/(A;) G whereas 
a H- or a W-solution is a distributional solution if v{k) G (see the Appendix I). 



In [7] the authors introduced a notion of solution that they call 'energy solution' 
(see the Appendix I). In fact an 'energy solution' is a W-solution which satisfies 
an additional property ensuring that = Wk (the additionnal property imposed 
is sufficient but not necessary to have this equality). Under this point of view an 
'energy solution' is slighty stronger than a W-solution. However, their existence 
result is obtained by assuming complicated conditions on the coefficients a and u 
which are not exactly satisfied in the physically relevant situation (Hp), but only in 
the following approximate situation: 



for some e > we have: 
a{s) = ai + a2\/s + e 



v{s) = ui + V2\Js + e 



On the contrary, our assumptions in Theorem 1 and Corollary 1 are very simple, 
and they are satisfied in (-f/p). 

Note also that we establish the regulartity property (2) which are not established in 
[7] (or in [10]). 

In the Appendix I we also give a new existence result for an 'energy solution'. 



In Theorem 2 we assume that (-ffo) and (-^2) hold. These assumptions are fulfilled 
in the physical situation [B.^ if ci2V\ = aiz/2. We then prove that u and k are Holder 



continuous. In particular we give here a positive answer to a central question put in 
[7] : k is bounded. Note that in this situation we clearly have Wk = Hk- 

We next establish the existence of a classical solution for Problem (P) by assuming 
some differentiability properties for a and v. These properties are fulfilled in the 
situation {H'^) if a2Vi = aiZ/2. 

It seems that this result is completely new: the existence of a classical solution for 
(P) was not studied in any previous work. 



1.4 Organization of the paper 



In the sequel n will always denote an arbitrary integer greater or equal to one, and 
C (possibly with subscript) will denote a positive real that does not depend on n, 
but that can differ from one part to another. 

We always consider the space Hq{Q) equipped with the gradient norm. 
The condition [Hq) is always assumed. 

• In section 2 we introduce an approximate sequence of solutions obtained 

by truncating the coefficients a and u. 

We immediatly obtain the basic estimates : 

/ z/„(A;„)|Vn„|'<C 
Jn 

Vp < - : ^ an{kn)Vkn ^ <C 
The point is that we establish the following fundamental estimates: 

||^n||L°°{n) < C 
J^an{kn)\Vkn\^<C {*) 



The first estimate above is proved by developping further a technique due to Stam- 
pacchia. 

The second is obtained under the assumption (Hi). The proof is based on the fol- 
lowing idea: if [u, k) is a solution of (P), we formally have ^ 

u{k)\Vu\'^ = - div {u{k)Vu).u+ div{u{k)uVu). (5) 

=fu 

In other words one can hope that the second member in the second equation in (P) 
is more regular than it seems. 

^ We thank Michel Chipot for this remark 



In fact, we prove that a similar relation to (5) holds for the approximate sequence. 
By using next that is uniformly bounded in L°°{Q), we obtain (*) which is the 
key estimate to prove Theorem 1. 

• In section 3 we extract from a subsequence converging to some element 

denoted by {u, k). Under the assumptions (-f^o) and [Hi), we directly obtain that 



• In section 4 we pass to the limit in the approximating Problems. In a first step 
we prove that m is a IV^-solution of (P).l. To do this, we use the test functions 
V = hg{kn)f (where ip G Wk H L°°{Q) and (hg) is a sequence of functions that cut 
off the large values), and we pass to the limits n — > oo, q oo. 

We next prove that the energies of the approximating sequence converge to the en- 
ergy ^z/(fc)|VMp. 

Finally we can pass to the limit in the second equation in order to prove that k is 
a distributional solution of (P).2. We then obtain Theorem 1 and Corollary 1 follows. 

• In section 5 we assume that (i^o) and {H2) hold. In a first step we obtain the 
estimate 



Hence k G L°° and by using the De Giorgi-Nash Theorem we prove the Holder con- 
tinuity of u and k. 

Next, by assuming additional regularity on u, dQ and / we can apply the Schauder's 
estimates and we prove Theorem 2. 

• In the Appendix I we study some relations between the notions of ly-solution, 
i?-solution, distributional solution, renormalized solution and 'energy solution' for 
Problem (P). We continue the discussion begun in Subsection 1.3 and we also es- 
tablish a new existence result for an 'energy solution' for Problem (P). 
In the Appendix II we recall some basic properties of Holder continuous functions. 

2 Approximating sequence and estimates 

We assume that [Hq) holds and we set 



u G H^{n) n k G H^in). 



We prove that moreover we have: 




||^n||L°°{Q) < C. 



z/„(s) =T„(z/(s)) 
a„(s) =T„(a(s)), 



(6) 
(7) 



where T„ is the truncated function defined by T„(t) = min(?7,, t). 
We consider the Problem of finding (ti„, kn) G (ifQ(n))^ such that 



For any n > 1, Problem (P„) is well posed because a„, Un G L°°(M) and a^^, u:^^ G 
L°°(M) by construction. 

It is proved in [7] that a solution kn) exists for any n > 1. Moreover, the following 
basic properties were established: 

fcn>0 (8) 

/ z/„(A;„)|Vn„|^<Ci (9) 

<\. jj\an{kn)Vkn\ <C2 (10) 

We now establish 

Lemma 3 The sequence Un is uniformly hounded in the L°°{Q)-norm, that is, 

\\Un\\L^(n)<Cs (11) 

Before proving this lemma we point out that the assumption / G with r > | 

made in (Hq) implies that 

f G W-^'f(n), with p = > 3. (12) 

3 — r 

This last property is easy to prove by using the Sobolev injection Theorem. 
Proof 

We will obtain the estimate (11) by using the technique presented on p. 108 in [12]. 
In order to prove that C3 is independent of n we have to detail the technique of 
Stampacchia. 
Let 

hn{u,v) := / Un{kn)VuVv. 



n 



Recall that / satisfies (12) and then by using a classical result (see [3]) there exists 
g G {LP{n)Y such that -dw{g) = f and \\g\\(LP(n)y^ < C\\f\\Lr{n). 
Hence the sequence u„ satisfies 



bn{un,v)= [ gWv yveH^in). (13) 
For s > 0, we define the measurable set A„(s) C f2 by setting 

An{s) = |x G : \Un{x)\ > s|. 



We also introduce 

(yj := max — s, 0) sgn{ur, 
It is proved in [12] that (f G Hq{Q) and 



(14) 



V(p = 'Vun in An{s) 
V(p = in fi\A„(s) 

By testing (13) with v = (p, we obtain 



A„(s) 



(15) 



Remark now that assumption z/(s) > 5 > in {Hq) imphes that Vnikn) > min((5, 1). 
Consequently the bilinear form 6„ is uniformly coercive on Hl{VL). By using this 
property together with the Holder inequality, we obtain from (15): 



< c 



An{s) 



I \A^'^\\ II 



Hence by using the Cauchy inequality together with the Holder inequality we obtain 



< Ci\\g 



|2 

\LP{n) 



p-2 



On the other hand, the Poincare-Sobolev inequality gives 



g\ 1/3 ~ 2 



(16) 



(17) 



Let now t > s. It is clear that An{t) C An{s) and consequently 



We set 



Anis] 



t - s\^]'^' > \t - S\' 



Vs > 



An{t) 



1/3 



For fixed n, ipn is a decreasing function, and from the estimates (16)-(18), we obtain 

Mt)<C3\Ms)fit-s)-' Vt>s>0, 

where we have used the notation P := ^^^^^^ > 1 and where C3 = C-s^Ci, C2, ||/||l'-)- 

Both quantity j3 and C3 do not depend on n. Hence by using Lemma 4.1 in [12] it 
follows: 

MS) = 0, 

1 /6 

where 6 = l^l^^ < 00 does not depend on n. 

This property tells precisely that (11) holds true with C3 = 9. □ 



Notice that the bilinear form 



(,u,v) 



an{kn)VuVv, 



is also uniformly coercive on Hq{Q). Moreover, the sequence 

hn := T„(z/„(fc„)|Vu„ 



is imbedded in L°°[Q). We can then apply again the technique of Stampacchia 
detailed in the proof of lemma 3, and obtain: 

forn > 1 : A;„ G (19) 

Nevertheless the control we have on is obtained from (9), which gives a uniform 
bound in the L^-norm for the sequence. This is not enough to obtain a uniform 
estimate for in the L°°-norm. 
However we can establish: 

Lemma 4 Assume that (Hq) and (Hi) hold. Then we have 

an{s) > 1ii^n{s), 7i = min(l, 7) (20) 

/ an{K)\Vkn\^<C5 (21) 
Jn 

Proof 

The estimate (20) is easy to obtain. Its verification is left to the reader. 

Let [un, kn) be the chosen approximating sequence. We have from (11) and (19) that 

Vn > 1 : M„, kn G H^{n) n L°°(fi) 

It follows (see [3]) that v := Un-kn G HI{VL) fl L°°(f2) is admissible for (P„).l and we 
get 2 

/ yn{kn)\VUn\^kn= / fUnkn- / h'n{kn)UnV U^V kn (22) 

JQ. Jn Jn 

By testing (Pn).2 with (p = kn, we obtain: 

/ an{kn)\Vkn\'^ = / T„ (/;;„) | Vm„ /i;„ < / z/„(A;„)|Vm„P/c„, (23) 
Jn Jn ^ ' Jn 

by using the properties T„(s) < s and (8). 
Hence, by combining (22) with (23) we have: 

/ := / an{kn)\Vknf < \ \fUnkn\+ Un{kn)UnVUnV kn (24) 

Jn Jn Jn, ^ ^ 

:=// :=III 

We can estimate the term II as follows: 



more generally: Un{kn)\Vun\'^ = -div(i/„(fcri)Viin).M„ +div(z^„(A;ri)ii„VM„) in 'D'{Q) 

^ V ' 



r Holder ineq. 

// < C3 / \fK\ < C4f\\LZ/2\\K\\L^ 

Jn 

Poincare-Sobolev Ineq. _ / f \ 1 /2 CJ^ / f \ 1 /2 

< C'lll/IL-v^l / |VA:„,P) < ^||/|L3/.( / a^{K)\VK\') 

•J ^ <J ^ 

Young Ineq. (J^ / 1 9 f r,\ 

< -^(-|I/IIl3/2 +e J^anikn)\Vkn\ ) for any e > given 

< \ [ a^iK)\WK\' + C2\\f Wls^. 

where 5 > is the constant given in (Hq). The last inequahty was obtained by 

~ ~ ~ 2 

choosing e = 5/(3Ci), using the estimate (11) and by setting C2 = 3Ci jS^. 



We next estimate the term III: 



111 = Un\l Un{kn)^Un\l l^niknj^K 



-1/2 



<o / «n(A;n)|V/i;„|2 + (74 / Vn{kn)\VUn\^, C^ = C^{C, 



3 Jn ' ' Jn 

where 6*3,71 are the constants that appear in (11) and (20). The last inequality 
follows from the Young inequality. 

Recall now the inequality (24) and use the estimates established for the terms II 
and III. We obtain: 



3 Jn 



- / an{kn)\Vkn\^ <C2 



L3/2(n) 



By using (25) together with (9) we finally obtain (21). 



(25) 
□ 



3 Basic convergence results for (n„, kn) 



The estimates established in the previous section allow us to extract a converging 
subsequence from (u„, /c„). We have 

Lemma 5 

1. Assume that (Hq) holds. Then we can extract a subsequence (still denoted by 
{un, kn) ) such that 



an{kn)V kn ^ a{k)V k in L^iVL), V < 7^ 



Un 
Un 



k a.e inVt 
u in Hq{Q) 
u inL'^iVL) 



(26) 

(27) 
(28) 
(29) 



2. Ij in addition the condition (H\) is fulfilled then we may assume that 

kn ^ k in Hq{Q) 



(30) 



Proof 

1. The properties (26) and (27) are obtained from (10). The property (28) is ob- 
tained by using the estimate (9) together with the assumption i/(s) > 5 > in {Hq). 
We estabhsh (29) from the estimate (11). 



2. By using Lemma 4 together with the assumption a{s) > 5 > in (Hq) we obtain 
(30). Notice that the k appearing in (26), (27) and (30) is necessarily the same in 
the three situations. □ 



We are able to prove additional regularity results for the element {u, k) introduced in 
Lemma 5. For technical reasons we introduce the sequence {/igjggM of real functions 
defined in [10] p. 185. It satisfies: 



\hg{s)\ < 1 W{q,s) e N X M (31) 
hg{s) = when \s\ > 2q (32) 

\h'(s)\ < - Vg e N, and a.e s e M (33) 
q 

hq ^ 1 uniformly on the compacts (34) 
Lemma 6 

1. Assume that (Hq) holds. Then the element {u,k) given in Lemma 5 satisfies 

I u{k)\Wu\^ < oo (35) 
Jo, 

2. Assume that in addition (Hi) holds. Then 

f a(k)\Vk\'^ < oo (36) 



Proof 

1. We take over the arguments presented in [10] p. 192. 
For g > 1, we set 

Vn,q ■= {hq{kn)Vn{kn)) ' Vm„ 

Let now q be fixed. The sequence \{hq{kn)i'n{kn)) ^ } >i uniformly bounded 
in L°°{Q). Consequently, {?7„,,g}n>i is bounded in {L?'{Vt)Y and we can extract a 
subsequence weakly convergent to some rjq G {L'^iVt))^. 
On the other hand, we have 



hq{kn)i'n{K)) ^ ^ {hq{k)vn{k)) ' a.e in Vt 



and thus rjg = {hq{k)v{k)^ ^ Vu. 
We now use a classical property of the weak convergence in LF'iVL): 



||'7j|L2{n) <liiminf ||?7„,g||L2(n) < li^inf { j^Vn{kn)\Vun\^) ' < C]'"^ 

where Ci is a constant independent of q given in (9). 
By using properties (34) and (31) we can see that 



ril z/(/i;)|VMp a.e. in ^2 
Vq< jy{k)\Vu\'^ 

Hence by the Fatou Lemma we finally obtain: 

J^u{k)\Vu\'^ < liminf llr/qll^a < Ci 

2. If the additional assumption (Hi) holds, then we have the estimate (21) and the 
previous reasoning allows us to obtain (36) □ 

4 The proof of theorem 1 



In the previous section we have proved that under (Hq) we can extract a converg- 
ing subsequence of {un,kn)- If moreover (Hi) holds then the limit {u,k) obtained 
satisfies: 



ueWkHL'^in) (37) 
k E H^in) (and in fact k E Wk) (38) 



4-1 Passing to the limit in (Pn)-l 

We recall that the space Wk was defined by 

Wk = \vE HliSl) : [v]k < oo 

We now establish: 



Lemma 7 Assume that (Hq) and (Hi) hold. Then the element {u,k) given in 
Lemma 5 satisfies (37), (38) and: 



[ v{k)VuVv = [ fv WveWk (39) 
Jn Jn 

Proof 

Let n > 1, g G N and if G Wk fl L°°{Q). We consider the function v := hq{kn)^p- By 
recalling the properties (31)-(34) of kg, we can verify that hq{kn) G HliVL) nL°°(f2). 
Consequently v G HliVt) fl By testing (P„).l with v, we obtain: 

/ := / Un{kn)hq{kn)VUnVip + / k' {kn)yn{kn)V UnV kn^ = / fhqkn^P (40) 

Jf7 



=// :=m 



In a first step we fix q and we study the behaviour of terms I, II and III when n 
tends to infinity. 

By using the property (32) we see that 

\i^n{kn)hq{kn)\ < uiax z/(s) := Cg, 

se[o,2g] 

and by using (32) together with (27) we obtain 

Vn{kn)hq{kn) — > v{k)hq{k) a.e in Vt. 

Consequently 

iyn{kn)hq{kn)Vip ^-^^u{k)hq{k)Vip in {L^{Q))^, 
and by also employing (28) we get: 

I ^ f iy{k)hq{k)VuV(p (41) 
We now estimate II. From (33) we obtain: 



//<! 

q J{q<k„<2q} 



h'n{kn)VUnVk 



< lblU--( / Vn{k^Wu^ff^\ f a„(fc„)|VA;„|2)'/' < -, (42) 

where the second inequality is obtained by using (20). 
For the last term we get 

/// ^ I fhq{k)ip (43) 



n^oo 



By using the estimates (41)-(43) together with (40) we obtain that for any fixed 
V? G Wfc n L°°(l]) the following holds 

/ u{k)hq{k)VuV^= I fhq{k)^+0{-). (44) 
Jn Jn q 



--Jl :=J2 



We next remark that the integrand in Ji converges for a.e. x E Qto u{k)'Vu'Vip when 
q tends to infinity. Moreover by using (31) together with the fact that (/? G Wk we can 
see that the integrand in Ji is dominated by \i'{k)'Vu'V(p\ G L^{Q). Consequently, 
by the Dominated Convergence Theorem we get 

Ji ^ I i^(kWuV(p 
9-^00 Jn 

Similarly we can see that 

9^00 Jn 

At this stage we have proved that 

/ u(kWuV(p = [ f(p ^ipeWkH L°°{n), (45) 
Jn Jn 

and it remains to show that the condition if G L°°{Q) is not necessary. 
Let ip G Wk and i eN. We consider ipi G Wk nL°°{fl) given by ipi = Ti{(p). By using 
some basic properties of Tj (see [7]), we see that \ipi\ < \ip\, \Vipi\ < \V(p\, fi ^ ip 
a.e , and Vipi — a.e in fl. Consequently, if we take ipi as test function in (45), 
we can pass to the limit ? — 00 and we obtain (39). □ 



In Lemma 7 we have showed that m is a VTfc-solution of (P).l. In order to prove 
Theorem 1 we have to prove that /c is a distributional solution of (P).2. We need 
first to establish: 

Lemma 8 Assume that (Hq) and (Hi) hold. Then, in addition to the results pre- 
sented in Lemma 5, we may assume: 

Un{kn)\^Un\'^ z/(fc)|Vu|2 mL^(l]) (46) 



Proof 

We test (Pri).l with the function u„,. By using (28) we obtain: 

/ i.„(A;„)|Vn„|2 ^ I fu= j y{k)\Vu\\ (47) 
Jn Jn Jn 

where the latter equality is obtained by testing (39) with u. 

We set Tjn '■= ^Ji^n{kn)Vun and ?7 := ^ u{k)Vu. The relation (47) tells us that 

hnWi^in) Mi^in) (48) 

We can next take over the arguments presented in [10] Lemma 5.3.4 in order to 
obtain: 

Tin ^ V in{L'{n)f (49) 

n^oo 

Finally properties (49) and (48) imply that the convergence is strong in (49), and 
(46) follows. □ 



4-2 The Proofs of Theorem 1 and Corollary 1 



Assume that {Hq) and (-f^i) hold. In Lemma 5 we have extracted a subsequence 
{un,kn) which converges in a certain sense to an element {u,k). This element has 
the properties (37)- (38). Next we have established (39). 
Let now ip G C^{Q). By using (26) we get: 

/ aJknWknVLp / a(k)VkVip (50) 

We next remark that the property (46) ensures that 

/ TnU{kn)\Wu^\^)ip ^ / z/(fc)|VM|V (51) 

Recall that the sequence {un,kn) satisfies {Pn)-2. Then relation (50) together with 
(51) allows to take the limit in (P„).2. We get: 

/ a{k)VkVip= [ z/(A;)|Vm|V G C,°°(fi) (52) 
Jn Jn 

Thus (P).2 is fulfilled in the distributional sense. 

At this point we have obtained (37), (38), (39) and (52). The proof of Theorem 1 is 
complete. 

Assume now that the condition (3) in Corollary 1 is fulfilled. By using (38) together 
with the Sobolev Injection Theorem we get k G L^(f2) and thus v{k) G L}{VL). Then 
we can conclude the proof of Corollary 1 by using Proposition 9 in the Appendix I: 
(m, k) is a distributional solution of (P). 

5 The proof of Theorem 2 

We assume in this section that {Hq) and {H2) hold. 

In this situation all the results presented in section 2 and section 3 are valid. For 
technical reasons we slightly modify the definition of a„ by setting 

an{s) := IVnis), (53) 

where 7 > is the constant appearing in (-^2) and Vn is defined as before. 
We will now consider Problems (P„) modified by the new definition (53) of a„. Nev- 
ertheless the modification is very sligth, and all the results presented in the previous 
section can be recovered easly. The verifications are left to the reader. 



We now prove that we have the new estimate: 



(54) 



In order to prove this result we set 

Xn ■■= K + (55) 

and we remark that (Pn)-2 leads to 

/ a„(A;„)VxnV(/9 = / fun^ Vc/? G HI{VL). 
Jn Jn 

Recall that a„(A;„) > 7min(l,5) > 0, an{kn) G L°°(f2) and note that the sequence 
fUn is uniformly bounded in U^Q) with r > 3/2. These properties are sufficient 
(see the proof of Lemma 3) to get the estimate 

||Xn||L°°(f7) < C, (56) 

where C does not depend on n. 

The estimate (54) is finally obtained by using Lemma 3 together with (56). 
Consequently, in addition to the properties in Lemma 5 we may assume that 

K^k in L°°(fi). (57) 

We will now prove that 

u,ke C°'"(n) for some a E (0, 1). (58) 

Let A := u{k). We have A, A-^ G L°°(fi) and 

/ AVmV0= / /0 V0G//o(fi). (59) 
Jq Jn 

Recall also that / have the property (12). Hence we can apply the De Giorgi-Nash 
Theorem (see for instance [5] Prop. 6 p. 683 or [8] Th. 8.22 and Th. 8.29). We obtain 
that u G C°'°i(n) for some ai G (0, 1). We next set x ■= k+{-f/2)u^. Then x e H^i^) 
and we have 

/ -VxV0= / fu<p V0GiJo'(^)- (60) 
Jn 7 Jn 

By using the fact that u G L°°{Q) in (60), we can again apply the De Giorgi-Nash 
Theorem to get x G C°'"2(f2) for some a2 G (0, 1). Hence also k is Holder continuous, 
and (58) follows. 

Let a G (0, 1) be a generic parameter that can differ from one part to another. We 

assume now that Oil is of class C^-", / g C°'"(n) and u G C^'"(M+). 

We will prove the second part of Theorem 2 by iterating the Schauder estimates. 

We have A = i/(/c) G C°'"(f2) (see the Appendix II) and then, by applying the 

Schauder estimate (see [4] Theorem 2.7 p. 154) on (59) we getM G C^'"(fi). Similarily, 

from equation (60) we obtain x ^ C^'"(f2) and thus k G C^'"(r2). 

Hence (see Appendix II) A G C^'°(r2). By iterating again the Schauder estimates (see 

now Theorem 2.8 p. 154 in [4]) we obtain that u and k are in C^'"(f2). 

Finally we see that [u, k) is a classical solution of (P). Theorem 2 is proven. 



Appendix I: Some relations between the notions of weak solution 

We give here some relations between the various notions of weak solution: W- 
solution, if-solution, distributional solution, renormalized solution, 'energy solu- 
tion'. 



Comparison with renormalized solution 



Proposition 9 

1. Any W- or H -solution (m, k) of Problem (P) that satisfies in addition k G Hq{Q), 
is also a renormalized solution. 

2. If i^{k) G L^{Q) then any W- or H -solution of Problem (P) is a also a distribu- 
tional solution of (P). 



1. Let {u, k) be a W^-solution of (P). Then the conditions (5.2.1)-(5.2.5) in [10] chap. 5 
are satisfied. We have to prove that (5.2.6) holds. 

Let h G C^(R) and (p G C^{Q) be arbitrarily chosen. We set v := h{k)(j). Then 
V G L°°{Q) and Vv = h{k)V(j) + h'{k)Vk(f). Let M < oo be such that the support of 
h being included in [— M, M]. We have 



Hence v G Wk- By testing (39) with v we obtain the relation (5. 2. 6). a in [10]. 

We remark that v is also admissible in (52). This allows us to obtain the condition 

(5.2.6).b in [10]. Consequently {u,k) is a renormalized solution of (P). 

If we consider a if-solution {u, k) of (P) we can take over the previous argument 

because the function v is now in Hj^. 

2. If z/(A;) G L^{Q) then we have C^(r2) ^ Wk- In consequence a Wk- or a 

i^fc-solution of (P).l is also a distributional solution of this equation. Hence (/i, k) is 
a distributional solution of (P). □ 



1. The first point in Proposition 9 tells that the notions of H- or VT-solution are 
stronger that the notion of renormalized solution. This fact is coherent with the 
second point established in Proposition 9: a H- or VT-solution is a distributional 
solution if ^{k) G L^{VL) whereas a renormalized solution is only a priori a distribu- 
tional solution if v{k) G L°°(fi) (see [10] p. 185). 

2. if we have k G Hl{VL) and if v satisfies the growth condition (3) then v{k) G L}{VL). 



We have: 



Proof 




Remarks 



Comparison with 'energy solution' 



We have seen that when u{k) G L^{Q) then any W- (or H-) solution is a distri- 
butional solution. Moreover the notion of ly-solution coincides with the notion of 
iZ-solution iff Wk = Hk (see [14]). 

Some sufficient conditions to have this last equality were established in [14] and in 
[7], but necessary and sufficient conditions are not known. 

Let us consider the following condition: 



It was schown in [7] that the first condition in (R) together with the property 
^-1 g /,°°(]R) (which is assumed in (Hq)) implies that Wk = H^. 

In [7] the authors introduced the notion of 'energy solution'. They impose {Hq) as 
the basic assumption. Then an 'energy solution' (m, k) for (P) is in fact a VT-solution 
which satisfies [R). This implies that Wk = Hk- The energy solution is also a H- 
solution, and moreover a distributional solution (because the first assumption in [R) 
implies that v{k) G L}{VL)). 

We see then that the notion of 'energy solution' (in the sense of [7]) has the advantage 



of unifying various notions by putting us in the situation where \Ji'{k) G H^{VL). 
The disadvantage is that we have to impose more complicated conditions on the 
coefficients a and z/, in order to obtain a solution. In particular in [7] Theorem 2.1, 
the authors prove the existence of an 'energy solution' under the assumptions (-f^o) 
and (i^s) (see below). 



This condition is not verified in the physical situation (-f/p), but only in the approx- 
imate situation {H'^). 



(R) 






V G Ci(R+) 

3 C > and 7 > 1/2 such that: 
\y\s)\<C VsG[0, 1] 

< C.s~^ Vs > 1. 



In Theorem 1 we obtain a VF-solution under much simpler conditions which are sat- 
isfied by iiip)- This solution is a distributional solution under an additionnal simple 



assumption (see Corollary 1) which is again satisfied in {Hp). 



Note also that in the first part of Theorem 2 we prove that under the assumptions 
(Hq) and {H2) (which are satisfied in (Hp) if ali^2 = ^22^1), the functions u and k are 
Holder continuous. In particular u{k) G which implies that Wk = Hk, and the 
notions of if-solution, ly-solution, distributional solution and renormalized solution 
coincide in this case. 



In order to conclude this Appendix we give a last existence result. Let {H4) be the 
following condition: 




u G C\R+) 

3C>0 s.t. ^ <c Vs G M. 

1/1 Q \ 



We have: 

Proposition 10 Assume that (Hq), (Hi) and (H4) hold. Then the W -solution given 
in Theorem 1 is an 'energy solution' (in the sense of [7]). 

Proof 

We have assumed that (-f^o); (Hi) hold and consequently all the results presented 
in the sections 2, 3 and 4 can be recovered. 



Let (m, k) be the 1^-solution given by Theorem 1. By using (2) we see that the 
second condition in (R) is satisfied. Nevertheless we cannot directly conclude that 
h'(k) G H^{^1), but we can obtain a new estimate for the approximating sequence 
{un, kn). More precisely, we have: 



\\^fM\\min)<C. (61) 

In fact, by using the property that /c„ G Hq{Q,) fl L°°{il) together with u G C^(M"'") 
we obtain i^(/c„) G H^{fl) fl L°°{fl), with Viy{kn) = i''{kn)Vkn- Recall now that 
^n{kn) = Tnivikn))- Heuce we have 

{kn}^ kn- 
it follows that: 



'Vjunikn) 



v'{kn)Vkn _ 
^{u(kn)<n} I — — - — ^{u„{k„)<n} 



jy'ikn) 



by (20) 

< C'l{i^„(fc„)<n} 



2^Un{kr, 

jy'ikn) 



2'\J l^n{kn)Q'n{kn 



O'nikn)^ kn 



^n{kn) 



O-nik-nj^kn — C 



y'{kn) 

V{kn) 



Clnikn)^ kn- 



Hence, by using (21) we obtain 

\\^\l^n{kn)\\L^{n) < C. 

Moreover \J'uJ^kn) = \J z^(O) on dVt and thus we obtain (63) by using a Poincare 
inequality. □ 

Remark 

The hypotheses made in Proposition 10 are verified under assumption {H'^. In the 
hypotheses, we require only very weak growth condition at infinity for u. For in- 
stance (contrarily to the result presented in [7]) the Proposition 10 works if we have: 

v{s) =1/1 + l/2e^l^ a{s) = ai + a2e^^^ /5i < ^2- 



Appendix II: Holder continuity and composition 

Let A C M"' and a e (0, 1). We recall that the space C'''°(A) of Holder continuous 
(with exponent a) functions on A is defined by: 

CO'"(A) = |/ : A^M s.t. VxoGA: snv\l^^hil^ < oo] 
^ x€A \X - xor 

More generally, for any integer k, the space C'^'"(A) is the space of those / G C^{A) 
whose kth derivative is in C°'"(A). 



A first elementary result tells that the product of two Holder continuous functions 
is an Holder continuous function. More precisely we have (see relation (4.7) in [8]): 

Lemma 11 Assume that /i,/2 G C°'"(A). Then /1./2 G C°'"(A) 

In Section 5 we used a function defined as a composition of two Holder continuous 
functions. We needed the following result: 

Lemma 12 Let Q be a compact in M.'^ and a G (0, 1). We consider the following 
three conditions: 



(A) XeC\R) and keC^'^iU) 

(B) X G C°'°(M) and k G C\n) 

(C) AgC^'°(M) and A; G C^'"(n) 

We have: 



1. Assume that (A) or (B) is satisfyed. Then X{k) G C°'"((]). 

2. Assume that (C) is satisfyed. Then X{k) G C"^'"(f2). 



Proof 

1. In this situation we clearly have X{k) G C°(n) and 

Ml := sup|A;(x)| < oo. (62) 

Let 

\\ik{x)) - \{k{Xo))\ 

l[X,Xo) .— . . 

\x — Xol 

We want to prove that 

sup I{x,Xo) < oo. (63) 

• Assume that (A) holds. Then in addition of (62) we have: 

|A(t)-A(to)| , \k{x)-k{xo)\ 
M2 := sup ^--^ < 00 and M3 := sup ^--^ < 00. 

t,ioG[-Mi,Mi] \t-to\ ^^^ggf^ |x-Xor 

Consequently: 

\x — Xol 

Hence (63) is satisfyed. 

• Assume now that (B) holds. Then in addition of (62) we have: 

\k(x)-k{xo)\ , |A(t)-A(to)| 
M4 := sup < 00 and M5 := sup ^-j^ ^ < 00. 

In this situation we can estimate /( follows: 

|A(A:(x))-A(fc(xo))| \k{x)-kixor < ^ 

^ ' °^ - \k{x)-kixor ■ k-^oh ' ■ 

Hence (63) is again satisfyed. 



2. Assume that (C) holds and let jj, := A(A;). Clearly /i G C^(n) and V/x = X'{k)Vk. 
We remark that A' G C'''"(M) and k G C^'"(f2). We can then apply the first point 
of this lemma to obtain: A'(fc) G C°'"(n). Moreover Vfc G (C°'°(n))'^. Hence the 
product A'(fc)V/c is Holder continuous. □ 



References 



[1] C. Bernardi, T.C. Rebollo, M.G. Marmol, R. Lewandowski, F. Murat: 'A model for 
two coupled turbulent fluids Part III: Numerical approximation by finites elements', 
Num. Math. 98, No 1,(2004), 33-66. 



[2] L. Boccardo, T. Gallouet: 'Non-linear elliptic or parabolic equations involving 
measure date', J. Func. Analysis 87, (1989), 149-169. 



[3] H. Brezis: 'Analyse fonctionnelle, theorie et applications', Masson, Paris Milan 
Barcelone (1993). 

[4] Y.-Z. Chen, L.-C. Wu: 'Second Order Elliptic Equations and Elliptic Systems', 
Translations of Mathematical Monographs, vol. 174, American Math. Society, 
Providence Rhode Island (1998). 

[5] R. Dautray, J-L. Lions: 'Analyse mathematique et calcul numerique pour les sciences 
et techniques, volume 2, I'operateur de Laplace', Masson, Paris New York, (1987). 

[6] P. Dreyfuss: 'Higher integrabiblity of the gradient in degenerate elliptic equations', 
accepted for publication in Potential Analysis. 

[7] T. Gallouet, J. Lederer, R. Lewandowski, F. Murat, L. Tartar: 'On a turbulent system 
with unbounded eddy viscosities'. Nonlinear Analysis 52 (2003), 1051-1068. 

[8] D. Gilbarg, N. S. Trudinger: 'Elliptic Partial Differential Equations of Second Order', 
Springer- Verlag, Berlin Heidelberg New York, (1977). 

[9] O. A. Ladyzenskaja, N. N. Ural'ceva: 'Equations aux Derivees Partielles de Type 
Elliptiques', Dunod, Paris, 1968. 

[10] R. Lewandowski: 'Analyse mathematique et oceanographie', Masson Paris, (1997). 

[11] J. Lederer, R. Lewandowski: 'A RANS 3D model with unbounded eddy viscosity', to 
be published in Ann. IHP, An. non Lin. (2006). 

[12] G. Stampacchia: 'Equations elliptiques du second ordre a coefficients discontinus', 
Les presses de I'universite de Montreal, (1966). 

[13] F. Serra Cassano: 'On the local boundedness of certain solutions for a class of 
degenerate elliptic equations', Bollettino f/.M./ (7)10-B (1996), 651-680. 

[14] V. V. Zhikov: 'Weighted Sobolev spaces', Shornik Mathematics 189(8) (1998), 1139- 
1170. 



